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Abstract. The systematics of different approximations within the self-energy-functional theory 
(SFT) is discussed for fermionic lattice models with local interactions. In the context of the SFT, an 
approximation is essentially given by specifying a reference system with the same interaction but 
a modified non-interacting part of the Hamiltonian which leads to a partial decoupling of degrees 
of freedom. The reference system defines a space of trial self-energies on which an optimization of 
the grand potential as a functional of the self-energy £2[E] is performed. As a stationary point is not 
a minimum in general and does not provide a bound for the exact grand potential, however, it is a 
priori unclear how to judge on the relative quality of two different approximations. By analyzing the 
Euler equation of the SFT variational principle, it is shown that a stationary point of the functional 
on a subspace given by a reference system composed of decoupled subsystems is also a stationary 
point in case of the coupled reference system. On this basis a strategy is suggested which generates 
a sequence of systematically improving approximations. The discussion is actually relevant for any 
variational approach that is not based on wave functions and the Rayleigh-Ritz principle. 

1. INTRODUCTION 

Lattice models of correlated electrons such as the single-band Hubbard model [1, 2, 3] 
represent one of the central issues in solid-state theory. One reason for this strong inter- 
est is that the Hubbard model is one of the simplest but non-trivial models that allow for 
a benchmarking of new theoretical concepts. In the recent years, dynamical cluster ap- 
proaches to the Hubbard model and its variants have become more and more popular [4]. 
Contrary to techniques that are based on the Ritz variational principle and on the opti- 
mization of wave functions, dynamical cluster concepts not only give information on the 
static thermodynamic properties of a system but also on the elementary single-particle 
excitations. These approaches can be divided into two groups: (i) cluster extensions of 
the dynamical mean-field theory (DMFT) [5, 6] and (ii) variational cluster extensions of 
the simple Hubbard-I approximation [1]. 

The DMFT can be understood as a mean-field theory which neglects spatial corre- 
lations but which fully takes into account temporal fluctuations. This is reflected in 
the DMFT approximation for the self-energy S,;j((o) which is local in the site indices 
E,j((u) = 5/jE((a) but shows up a non-trivial and in general strong dependence on the 
excitation energy O). Spatial correlations are systematically restored in the dynamical 
cluster approximation (DCA) [7, 4] or in the ceUular DMFT (C-DMFT) [8, 9]. Contrary 
to the original DMFT, where the self-energy is generated from a model where a single 
correlated site is embedded into a continuous non-interacting medium ("bath"), the clus- 
ter extensions employ more complicated reference systems where the single-site impu- 
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FIGURE 1. Illustration of two reference systems H' generating two different approximations: The 
variational cluster approximation (VCA) and the cellular dynamical mean-field theory (C-DMFT). The 
reference systems consist of Lc > 1 correlated sites (non-zero Hubbard-f/) per cluster. The intercluster 
hopping (dashed lines) is switched off. In case of the C-DMFT an uncorrelated, continuous bath is attached 
to each of the original sites in addition. For both, the VCA and the C-DMFT, all one-particle parameters 
of H' are considered as variational parameters. 



rity is replaced by a cluster of Lc > 1 correlated sites. This generates self -energies with 
off-site elements. The bath parameters are determined by a so-called self-consistency 
equation which relates the reference (impurity/cluster) model with the original (Hub- 
bard) model. 

The Hubbard-I approximation represents a very simple approximation scheme which 
originally was constructed [1] by a more or less ad hoc decoupling of the equations of 
motion for the one-particle Green's function. Equivalently, however, it can be understood 
as a scheme which approximates the self-energy of the original Hubbard model by the 
self-energy of an atomic model consisting of a single correlated site (Lg — 1). From this 
perspective, a cluster generalization is straightforward and yields the cluster-perturbation 
theory (CPT) [10, 11]. The CPT can also be considered as the first non-trivial order in 
a systematic expansion in powers of the inter-cluster hopping parameters. The usual 
CPT uses a cluster of finite size > 1 which is cut out of the original lattice to 
generate the approximate self -energy. The main idea of the variational CPT (V-CPT) 
[12, 13] is to optimize the self-energy by varying the parameters of the cluster. This is 
reminiscent of the optimization of the self-energy by varying the bath parameters in the 
cluster extensions of the DMFT For the construction of a thermodynamically consistent 
approximation the variational aspect is essential [14]. It is therefore reasonable to call 
this a variational cluster approach (VCA). 

Both types of approximations, (i) and (ii), can be obtained within a general framework 
which is known as the self-energy-functional theory (SFT) [15, 16]. The SFT is based 
on a variational principle dQ.t^u[^] = [17] which goes back to the original ideas of 
Luttinger, Ward, Baym and Kadanoff in the sixties [18, 19] and which provides a very 
general framework to construct dynamical approximations. Let the original (Hubbard- 
type) model with one-particle and interaction parameters t and U be given on a lattice 
consisting of L sites (with L — > oo). Consider then a partitioning of the lattice into 
disconnected clusters with a finite number Lc of correlated sites (and possibly also 
a number of Lb additional uncorrelated bath sites attached to each of the correlated 
sites). The model on the truncated lattice (the "reference system") is therefore given 



by modified one-particle parameters t' and serves to define trial self-energies ^t',U 
for the variational principle. The trial self-energy is varied by varying the one-particle 
parameters of the reference system In this way one can search for a stationary point 
of the self-energy functional on the restricted space of self-energies defined by a simpler 
reference system: 

An,,j,[E,,j,]=0. (1) 

The type of the approximation is determined by the choice of the reference system, i.e. 
by the cluster size Lc, and by the number of bath degrees of freedom L),. The DMFT 
is obtained for Lc = 1 and Lb = for the C-DMFT one needs Lc > 1, and the VGA 
is specified by the choice Lg > 1 and Lb = (see Fig. 1). Clearly, there are more pos- 
sibilities. Approximations constructed in this way are dynamic and thermodynamically 
consistent in general: Via the self-energy at the stationary point, they provide informa- 
tion on the one-particle excitations and an explicit, though approximate, expression for 
a thermodynamic potential from which all static quantities of interest can be derived. 

Contrary to the Ritz variational approach, the SFT cannot predict exact upper bounds 
for the grand potential. From the Ritz principle, or from its generalization for arbitrary 
temperatures [20], one has Q.t u [p] > ^t.u^ i-^- the grand potential at an arbitrary density 
matrix p represents an upper bound of the exact grand potential O^t^jj. On the other hand, 
nothing prevents that Q-t^u < ^t,U for some t' within the SFT. 

This raises a number of questions which are addressed in the present paper: (i) Is 
there more than a single stationary point of the self-energy functional, i.e. is there more 
than a single solution of the Euler equation (1)? (ii) If this is the case, which one is to 
be preferred? (iii) Comparing two different approximations resulting from two different 
choices of the reference system, which one is more reliable? 

These are questions that refer quite generally to any variational principle that does 
not share with the Ritz principle the "upper-bound property". It will be argued that 
always taking the stationary point with the lowest SFT grand potential is a strategy that 
is generally unacceptable. A different strategy is suggested instead. 

The paper is organized as follows: The next section briefly reviews the basic concepts 
of the SFT. An explicit form of the Euler equation (1) is derived in section 3. Section 4 
presents an analysis of the Euler equation for the case of a reference system composed 
of two decoupled subsystems. This forms the basis for the general discussion on the 
relative quality of different approximations and a systematic way to approach the exact 
solution in section 5. The main conclusions are summarized in section 6. 



2. SELF-ENERGY-FUNCTIONAL THEORY 

The central idea of the self-energy-functional theory (SFT) is to make use of the univer- 
sality of the Luttinger-Ward functional 4>t/[G] [18] or of its Legendre transform Li/[E]: 
For a system with Hamiltonian H = Ho{t) +Hi (U), where t are the one-particle and U 
the interaction parameters, the functional dependence Fu[- ■ ■] is independent of t [17]. 
The grand potential of the system at temperature T and chemical potential jU can be 



written as a functional of E: 



at^uin = Trln(G^; -E)-i +Fu[n , (2) 

where Gf^ = {i(On + jU — is the free Green's function and Tr = ri^„e'^"*^ tr with 
the usual trace tr and the Matsubara frequencies a)„ = {In + 1 ) Trr for n = 0, ± 1 , ±2, 
After Legendre transformation, the basic property of the Luttinger-Ward functional reads 
as r^^5Fi/[Z]/5E = — G{/[E]. This implies that Dyson's equation can be derived by 
functional differentiation, 5Q.t,u\P\/ = 0- Hence, at the physical self-energy L = 
T.t u, the grand potential is stationary: 5Q.t^u[Lt,u] = 0- 
Due to the universality of /^t/P]' one has 

^t',u[n = Fuin +Tr\n{G;,yL)-' (3) 

for the self-energy functional of a so-called "reference system" which is given by a 
Hamiltonian with the same interaction part U but modified one-particle parameters t': 
H' = Ho{t') +Hi{U). The reference system has different microscopic parameters but is 
taken to be in the same macroscopic state, i.e. at the same temperature T and the same 
chemical potential ji.By a proper choice of its one-particle part, the problem posed by 
the reference system H' can be much simpler than the original problem posed by H, 
such that the self -energy of the reference system ^.f/ y can be computed exactly within 
a certain subspace of parameters t'. Combining Eqs. (2) and (3), one can eliminate the 
functional Fu [L] . Inserting as a trial self-energy the self -energy of the reference system 
then yields: 

^t,u[^t',u] = ^t',u + TMG;,o -^t',u)~^ -TrlnG,/,!/ , (4) 

where ^t',u and Gt'_u are the grand potential and the Green's function of the reference 
system. This shows that the self-energy functional can be evaluated exactly on the 
subspace of trial self-energies that are generated by the reference system. Solutions of 
Eq. (1) represent stationary points of the functional on this subspace. For further details 
of the approach and for different applications see Refs. [12, 13, 14, 15, 16, 17, 21, 22, 
23, 24, 25, 26, 27, 28, 29, 30]. 



3. SFT EULER EQUATION 

Let {\oc)} denote the orthonormal set of one-particle basis states. Then tajj are the 
elements of t, Gf U-apiicOn) — ((ca;c^))r,i/ are the elements of Gt^u = Gt,u{io^)-, etc. 
Carrying out the partial differentiation d / dt' in Eq. (1), one arrives at 

f ^ rr.\\ ^^t',U;al5{ia>n) 

Note that there are as much (non-linear) equations as there are unknowns t'^^ . The Euler 
equation (5) can be derived from the representation (4) for the SFT grand potential 



by using the (exact) relation dQ-tt u/dt'^^ {c^aCp)t',U = TY,,^e'°^"^^ {icOn + - t' — 

'Lff jj)p^. The Euler equation is trivially fulfilled for t' = t since (G^q — i/)"^ = Gf ij. 

In all practical situations, however, the point t' = t does not belong to the parameter 
space of the reference system since the t' must be chosen such that the problem posed 
by H' is exactly solvable. 

The term d'Lfi u.fjp{iC0n)/dt' may be considered as a projector. In the space of self- 
energies, <5Ef',i/;aj8 (io^n) I dt' is a vector tangential to the hypersurface of t' representable 
trial self-energies Et',f/- Hence, the Euler equation Eq. (5) determines the self-energy 
from its exact conditional equation (Dyson's equation) but projected onto that hypersur- 
face by taking the scalar product with the projectors dTLti u / dt'. 

The projectors can be determined more explicitly by carrying out the t' differentiation. 
Writing E = E,/ f/ and G' — Gf'^u foi" short, one has 

^ap (iO^n) = iia)n + IJ')Sap - t'ap - G'~p {i(On) (6) 

from Dyson's equation for the reference system. Hence: 

—d7' --^ocp5p^ • (7) 

Making use of the relation 5{B- A)7.^/5Amn = (B - (B - A)^^^ which holds for 
two not necessarily commuting matrices A and B, 

■£} = -0ap0iio+2^G a^{l(0„) g- Gyp{lCO„). (8) 

"'per IJ.V "'■pa 

To calculate the linear response of the Green function when varying the one-particle 
parameters, the S matrix shall be introduced by the definition 

St:,u{\/T) = cxpiHiiU)/T)expi-iHoit')+HiiU)-iiN)/T) . (9) 

is the particle number operator. Note that as compared to the conventional definition 
for S, the roles of the free and of the interaction part of the Hamiltonian are interchanged, 
i.e. HQ{t') — ilN is considered as a perturbation here. With the time ordering operator 
the S matrix can be written as 

5,.t/(l/r) = ^exp|^-^'^^tfi:pr;^-AiV)^(^)^<T(^)j • (10) 

Here the notation 0{t) = exp(//i(f/)T)Oexp(— //i(t/)T) is used: The imaginary time 
dependence is due to Hi {U ) only. 

Now the Green's function G^y(/(0„) — Jq <iTe'""'^G|^y (t) can be written as 

tr.-^i(^)/^^5c^(T)c|;(0) 



FIGURE 2. Diagrammatic representation of the relation between the two-particle Green's function L' 
and the two-particle self-energy T' . Double lines represent the one-particle Green's function G' . 



with S = Sfi u{\/T) for short. Again, the (interchanged) interaction representation is 
used with the t dependence of c^{x) being due to H\{U). The t' dependence of the 
Green's function is due to S only: 



7 S /■ i / i 

— = -^/ JtSJ(t)?^(t)5 (12) 

Using this in Eq. (11), 
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pa -^0 



(13) 



where L' is a two-particle dynamical correlation function of the reference system H': 

4iP2P3P4(^1''^2;T3,T4) = {^Cp,{Ti)cp^{T2)cl^{T3)cl^{r4)) 

-(^Cp,(Ti)4(T4))(^Cp,(T2)cj3(T3)) (14) 

and T+ = T + 0+. Here the average and the time dependence is due to the full Hamilto- 
nian: 0(t) = exp({H' - fiN)T;)Oexp{-(H' - fiN)^). Defining 



one has 



and thus 



— = -^MC7pv(«6^) (16) 



-pa 



, - -^ap^fia-Y^G' a^{(0)lJ^^p^{(0)G\p{(0) . (17) 



Introducing the two-particle self -energy of H' (see Fig. 2), 



^P.P.PspM^^) = f 5G'p^p^{a>2) ' ^^^^ 



where the functional of G' is the functional derivative of the Luttinger-Ward functional, 
this can also be written as 

= -T'L^LKy^piifOr^ML'^apviiOirn) , (19) 
^^pa m nv 




FIGURE 3. Diagrammatic representation of the Euler equation (20). All variables except for p and a 
are summed over. Rectangular box: two-particle self-energy F'. Circle and adjacent double lines: two- 
particle response function L'. Double line: Green's function of the reference system G'. The solid line 
denotes the approximate Green's function (G^q — E)~' of the original system. 



which finally yields the Euler equation in the form (see Fig. 3): 



Gf^{iC0n) -T.{ico„) 



- G'{ia)n) ^'av^lp (^"^' i(^)L'^apv • 
- jSa 

(20) 

If the system H' not only consists of one-particle orbitals belonging to H but also 
includes additional (uncorrelated) bath orbitals, one has to be careful with the orbital 
indices. Throughout the above derivation, a,j8,7, 5 refer to the orbitals of H, while 
li,v,p,G refer to the orbitals of H'. An orbital index ju of H' runs over n = a = 
1,2, ...,M (the correlated orbitals of H' which are identified with corresponding orbitals 
of H, where usually M — oo) and additionally over /i=M+l,M + 2,... (the uncorrelated 
bath orbitals). G' denotes the Green function of the reference system with the elements 
G'^y{i(On) ■ On the correlated orbitals n = a, v = j5 one has G'^^ (iWn) = Gu-^^^ {i(On) P] • 
Recall that Gu [£] is the inverse functional of Ei;[G] which only includes the propagators 
between correlated sites a and j8 . When additional uncorrelated sites are considered, 
the equation (6) is not the complete Dyson equation in H' but only the block with a, j8 

elements. Note that G' ^ means matrix inversion with respect to all orbitals of H' . 

As compared with the DMFT or the C-DMFT self -consistency equation, the SFT Eu- 
ler equation (20) is more complicated as it involves dynamical two-particle correlation 
functions of the reference system. As the reference system is assumed to be exactly solv- 
able these are accessible, in principle. For practical calculations, a modified version of 
the Euler equation has been suggested [22] and shown to allow for an extremely precise 
determination of a stationary point of the self-energy functional. For the purpose of a 
general discussion, however, the form (20) is more useful. 



4. A THEOREM ON DECOUPLED REFERENCE SYSTEMS 

Let a reference system consist of two subsystems A and B. Subsystem A is defined as 
the set of orbitals |oca)> and subsystem B is given by the rest of the orbitals la^), i.e. 
the complete (and orthonormal) one-particle basis is (Icca), |o;b)}. Typically, A and B 
are given by two disjoint sets of sites. The Hamiltonian of the reference system can be 
written as 

H'^H'^+H's+H'j,s, (21) 



where //^ only acts in the Fock space of {Icca)} and H'g in the Fock space of 
Hence, the commutator [//^,//^]_ = 0. //^^ is a term which couples the dynamics of the 
two subsystems and is assumed to be a one-particle operator, i.e. a coupling due to a two- 
particle interaction part of the Hamiltonian is excluded. This is satisfied, for example, 
in case of the Hubbard model if the subsystems are given on disjoint sets of sites as 
the Hubbard interaction is local. The coupling term can then be written as //^^ = //o(^) 
where V is the matrix of one-particle coupling parameters Va^ae with Ua&A and as e B. 

A given reference system specifies a certain space of trial self-energies for the self- 
energy functional and thereby a certain approximation. What is the relation between an 
approximation given by the reference system//' — H^+Hg +H^b approximation 
given by the decoupled system H" = H'^+ H'gl With the above preconditions, the 
following theorem holds: Any stationary point of the self-energy functional on the 
subspace of self-energies defined by the decoupled system H" is also a stationary point 
on the subspace of self-energies defined by the coupled system H', namely at V = 0. 
Writing H' = Ho{t')+Hi{U) = Ho{t") + Ho{V) + Hi{U) and H" = Ho{t") +Hi{U) 
and 



the theorem is: 



d d 



= 0. (23) 



While the theorem is not trivial, it complies with intuitive expectations: Going from 
a more simple reference system to a more complicated reference system with more 
degrees of freedom coupled, should generate a new stationary point with V ^ 0; the 
"old" stationary point with V = 0, however, is still a stationary point in the "new" 
reference system. Therefore, coupling more and more degrees of freedom, introduces 
more and more stationary points of the self -energy functional, and none of the old ones 
is "lost". 

For the proof the results of the preceeding section are needed, in particular the 
representation of the projector d'Lccji{iO)n) / dt'p^ in Eq. (17). One has to distinguish 
between the following different cases: 

(i) a,j8 e A and p,a e B: For V = the Green's function as well as its inverse does 
not couple orbitals of different subsystems, e.g. G'^j^{iC0n) = if a E A and /i G 5. 
Hence, in Eq. (17) there can be non-zero contributions for /i, v G A only. Since p, C7 G 5 
and V = 0, the first term of the two-particle Green's function in Eq. (14) decouples and 
thus: 

- (^C^(T)4(0))(^C,(T')ct«)) 
= (^C^(T)4(0))(^Ca(T')cJ«)) 

- (^c^(T)4(0))(^Cc,(i:')4«))=0' (24) 

which implies d'L(^p{iG)n) / dt'p^ = 0. The case E B and p,<7 G A can be treated 
analogously. 



(ii) (X,j5,p G A and o E B:ln Eq. (17) there can be non-zero contributions for /i, v G A 
only. For V — the Green's function L' decouples and vanishes since (C(j(t')) = for 
fermions. Consequently, d'L(^p{i(On) / dtp^y = 0. The same type of reasoning applies to 
the cases a,P,o eA, p e B and a,o,p eA, P e B and j8,0',pGA, aG5as well as to 
those cases with the roles of A and B interchanged. 

(iii) a,p G A and j8, a G fi: In this case, ;U G A and v eB which implies (for V = 0) 
that the second term of the two-particle Green's function in Eq. (14) vanishes and the 
first term decouples: 

L'^opyi^.^'-^^.O) = (^c^(t)c^(tOcJ(t;)4(0)) 

= -(^c^(t)c;(t;))(^c^(04(o)) 

= -G'^p{'^-T')G'„y{T'). (25) 

This yields: 

L'^^opviicon) = - fj^ fj^ dxdx'e'"^'G'^p{x-x')G'^^{x') 

= -G'^p{ia)n)G'^^{ia)n) (26) 

and thus 

^^fj'^^ = -5ap5^^ + £G'„i(/(0„)G;,p(;(0„)G'^^(/«„)G';^(/w„) =0. (27) 

Analogously, the projector vanishes if a, p G 5 and /3, C7 G A. 

(iv) In the case a, cr G A and I5,p eB and, analogously, for a, cr G 5 and /3 , p G A, one 
is led to anomalous correlation functions of the form (cc) which vanish if spontaneous 
U(l) symmetry breaking is excluded as it is done in the derivation of the Euler equation 
(1) from the very beginning. As a consequence, one has d'L(^^{iGyn) / dt'p^ = in this 
case, too. 

Thereby all possibilities have been enumerated with the exception of the two 
cases a,/3,p,a G A and a,/3,p,(7 G B. Here, there is no reason for the projector 
dlLa^ (icOn) I dt'pfj to vanish even if V = 0. 

These last two cases in fact correspond to variations on the space of trial self-energies 
given by the decoupled system H". If there is a stationary point ^{t") on this smaller 
space, this must necessarily represent a stationary point on the larger space given by 
H', too: Namely, in the additional cases to be considered, the Euler equation is fulfilled 
trivially since, as shown above, the projector d'Lo,p{ico„) /dtp^y vanishes. Summing up, 
this shows that any stationary point of the self-energy functional on the smaller subspace 
is also a stationary point on the larger subspace of the coupled reference system, namely 
at V = 0. This proves the theorem. 

Put in another way, the theorem states that as a function of a parameter (set of 
parameters) V coupling two separate subsystems, 

^tM^"+v,u] = ^t,u[^,"+o,u] + eiy^) , (28) 

provided that the functional is stationary at 'Lt'^u when varying t" only (this restriction 
makes the theorem non-trivial). 
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FIGURE 4. Different possible reference systems with the same interaction as the single-band Hubbard 
model on a square lattice. Filled circles: correlated sites with U as in the Hubbard model. Open circles: 
uncorrected "bath" sites with f/ = 0. Lines: nearest-neighbor hopping. Big circles: continuous bath 
consisting of Lb = °° bath sites. Reference systems H'^, H'^, H'^ generate variational cluster approximations 
(VCA), //g a dynamical impurity approximation, H[ the DMFT, and H'^ an intermediate approximation 
(VCA with one additional bath site per cluster). 



5. HIERARCHY OF STATIONARY POINTS 



To be explicit and to simplify the discussion, the single-band Hubbard model 

^ = -^ E 4c;cT + t/£n,rn,i (29) 

on a two-dimensional square lattice with nearest-neighbor hopping is considered in the 
following. Nevertheless, the discussion is completely general and applies to arbitrary 
correlated fermionic lattice models with local interactions. A possible reference system 
U' must have the same local Hubbard interaction; the hopping part, however, can be 
modified arbitrarily. A number of different reference systems are shown in Fig. 4. 

To discuss a first consequence of the theorem, one should distinguish between "trivial" 
and "non-trivial" stationary points for a given reference system. A stationary point is 
referred to as "trivial" if the one-particle parameters are such that the reference system 
decouples into smaller subsystems. If, at a stationary point, all degrees of freedom (sites) 
are still coupled to each other, the stationary point is called "non-trivial". In fact, all 
numerical results that have been obtained so far [12, 13, 14, 15, 16, 17, 21, 22, 23, 24, 
25, 26, 27, 28, 29, 30] show that there is at least a single non-trivial stationary point for 
any reference system. 

Once this is assumed to be true, then, for a given reference system, there must be 
several stationary points. Consider the reference system H'^ in Fig. 4, for example. Here, 



there are four intra-cluster nearest-neighbor hopping parameters which are treated as 
independent variational parameters. A non-trivial stationary point would be a stationary 
point with t[ , ?2, ^3, ^4 7^ (or , ^2' ^3 7^ ^' ^4 = ^ second stationary point is then found 
for = — and some , 7^ since, according to the theorem, this corresponds to 
a non-trivial stationary point generated by the reference system Another stationary 
point is obtained with t[ — — t'^ — — since this corresponds to a stationary point 
generated by H!^. (Note that the one-particle energies e- = t'- are variational parameters, 
too.) This shows that within a given approximation, i.e. for a given reference system, a 
non-trivial stationary point has always to be compared with several (on that level) trivial 
stationary points. 

Now, it is important to note that a stationary point of the self-energy functional 
£lf u\Lti u] is not necessarily a minimum. In general, a saddle point is found. This 
is demonstrated, for example, by the calculation in Ref. [15]. Furthermore, there is 
no reason why, for a given reference system, the SFT grand potential at a non-trivial 
stationary point should be lower than the SFT grand potential at a trivial one. And 
finally, it cannot be ensured that the SFT grand potential, evaluated at a given trial self- 
energy, is always higher than the exact grand potential, i.e. fl^ ^/ [£;/ (/] < Q.t u may 
be possible. This stands in sharp contrast to the Ritz variational principle. The fact 
that the spectrum of the Hamiltonian (after a constant energy shift) is always positive 
definite guarantees the upper-bound property (^|//|^) >Eq. It is not surprising that 
this upper-bound property is lost within the SFT as the approach does not refer to 
wave functions |^) at all. This is probably characteristic for any dynamical variational 
approach, i.e. for variational approaches based on time-dependent correlation functions. 
Green's functions, self-energies, etc. 

Consider the case where there is a non-trivial stationary point and a number of trivial 
stationary points for a given reference system. Despite the above reasoning, an intuitive 
strategy to decide between two stationary points is to always take the one with the lower 
grand potential Q.t;u^t',u\- A sequence of reference systems (e.g. H'^, H^, H'^, ...) in 
which more and more degrees of freedom are coupled and which eventually recovers the 
original system H itself, shall be called a "systematic" sequence of reference systems. 
For such a systematic sequence, the suggested strategy will produce a series of stationary 
points with monotonously decreasing grand potential. This is reminiscent of the Ritz 
principle. Furthermore, by comparing the trends of the SFT grand potential for two 
stationary points as functions of an external parameter, one can easily identify level 
crossings as well as continuous or discontinuous phase transitions and interprete them 
consistently within the framework of equilibrium thermodynamics. 

Unfortunately, however, the strategy is useless because it cannot ensure that a system- 
atic sequence of reference systems generates a systematic sequence of approximations 
as well: Within the SFT, one cannot ensure that the respective lowest grand potential in a 
systematic sequence of reference systems and corresponding stationary points converges 
to the exact grand potential. This means that despite the fact that the complexity of the 
reference systems increases, the stationary point with the lowest SFT grand potential 
could be a trivial stationary point, i.e. could be associated with a very simple reference 
system only (like H'^ or for example). Such an approximation must be considered 
as poor since the exact conditional equation for the self-energy is projected onto a very 
low-dimensional space only. 



a) b) c) d) 




FIGURE 5. Schematic trends of the SFT grand potential i2 as a function of a variational parameter V 
which is assumed to couple two different subsystems of the reference system, i.e. V = corresponds to 
the decoupled case and (according to the theorem of Sec. 4) must always represent a ("trivial") stationary 
point. Circles show the stationary points to be considered. According to Rl, the trivial stationary point 
V —0 has to be disregarded in all cases except for a). The arrow marks the respective optimum stationary 
point. Note that in case c) the SFT grand potential is higher for the non-trivial as compared to the trivial 
stationary point. In case d) the rule R2 applies. 



Therefore, one has to construct a different strategy which necessarily approaches the 
exact solution when following up a systematic sequence of reference systems. Clearly, 
this can only be achieved if the following rule is obeyed: 

• A non-trivial stationary point is always preferred as compared to a trivial one (RO). 

A non-trivial stationary point at a certain level of approximation, i.e. for a given refer- 
ence system becomes a trivial stationary point on the next level, i.e. in the context of a 
"new" reference system that couples at least two different units of the "old" reference 
system. Hence, by construction, the rule RO implies that the exact result is approached 
for a systematic series of reference systems. 

Following the rule (RO), however, may lead to inconsistent thermodynamic interpre- 
tations for the case that a trivial stationary point has a lower grand potential as the non- 
trivial one. To avoid this, another rule is necessary: 

• Trivial stationary points have to be disregarded completely unless there is no non- 
trivial one (Rl). 

This automatically ensures that there is at least one stationary point for any reference 
system, i.e. there is at least one solution at any level of the approximation. Clearly, Rl 
makes RO superfluous. 

To maintain a thermodynamically consistent picture in case that there are more than a 
single non-trivial stationary points, one needs the following rule: 

• Among two non-trivial stationary points for the same reference system, the one 
with lower grand potential has to be preferred (R2). 

The rules are illustrated by Fig. 5 which gives different examples. Note that the grand 
potential away from a stationary point does not have a direct physical interpretation. 
Hence, there is no reason to interprete the solution corresponding to the maximum in 
Fig. 5, c) as "locally unstable". The results of Ref. [21] (see Figs. 2 and 4 therein) nicely 
demonstrate that with the suggested strategy (Rl, R2) one can consistently describe 
continuous as well as discontinuous phase transitions. 



It should be stressed that the above rules Rl and R2 are unambiguously prescribed 
by the general demands for the possibility of systematic improvement and for thermo- 
dynamic consistency. There is no acceptable alternative to this strategy. Note that the 
strategy reduces to the standard strategy (always taking the solution with lowest grand 
potential or, for T = 0, the lowest ground- state energy) in case of the Ritz variational 
principle because here a non-trivial stationary point does always have a lower grand 
potential as compared to a trivial one. 

There are also some consequences of the strategy which might be considered as 
disadvantageous but must be tolerated: (i) For a sequence of stationary points that 
are determined by Rl and R2 from a systematic sequence of reference systems, the 
corresponding sequence of SFT grand potentials necessarily converges to the exact 
grand potential but not necessarily in a monotonous way. For example, the exact grand 
potential might be approached from below or in an oscillatory way. (ii) Given two 
different approximations specified by two different reference systems, it is not possible 
to decide which one should be regarded as superior unless both reference systems belong 
to the same systematic sequence of reference systems. In Fig. 4, for example, one has 
H'^ < H^^ < < H'^ where "<" stands for "is inferior compared to". Furthermore, 
< and < but there is no relation between //^ and H^, for example. 

6. CONCLUSIONS 

A dynamical variational principle is a principle of the form 5^2 [X] = where ^2 is a 
thermodynamic potential and X is a dynamical quantity that refers to excitations of the 
system out of equilibrium but in the linear-response regime. A common characteristic 
of the different dynamical variational principles used in solid-state theory [16] is that 
stationary points are saddle points rather than minima in general and that the thermody- 
namic potential at a stationary point cannot serve as an upper bound of the true potential. 
One of the most famous approximations that can be constructed in this context is the dy- 
namical mean-field theory and, in fact, there is no general proof (for finite-dimensional 
lattice models) that ^^dmft > ^exact so far. 

Having these problems in mind, it becomes questionable how to judge on the relative 
quality of two different approximations resulting from two different stationary points 
of a dynamical variational principle. It has been shown here that at least within the 
context of the self-energy-functional theory there is an answer to this question which is 
prescribed by demanding approximations to be thermodynamically consistent as well as 
systematic and which is summarized by the rules Rl and R2 in Sec. 5. It has turned out 
that the intuitive strategy of always preferring the stationary point with the lowest SFT 
grand potential is unsystematic and therefore unacceptable. The essence of the correct 
strategy, on the other hand, is to disregard, as far as possible, those stationary points that 
(at a certain level of the approximation) are trivially induced due to a partitioning of the 
reference system into subsystems with fully decoupled degrees of freedom. 
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